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Objectives 

This recently completed three-year project was undertaken by the PI at the University of Florida, 
NASA Goddard and JPL, and by the Co-I and Collaborator at the University of Florida. The 
funding was used to support a continuation of research conducted at the University of Florida over 
the last decade which focuses on the dynamics of dust particles in the interplanetary environment. 
The main objectives of this proposal were: 

• To produce improved dynamical models of the zodiacal cloud by performing numerical sim- 
ulations of the orbital evolution of asteroidal and cometary dust particles. 

• To provide visualizations of the results using our visualization software package, SIMUL, 
simulating the viewing geometries of IRAS and COBE and comparing the model results with 
archived data. 

• To use the results to provide a more accurate model of the brightness distribution of the 
zodiacal cloud than existing empirical models. 

In addition, our dynamical approach can provide insight into fundamental properties of the 
cloud, including but not limited to the total mass and surface area of dust, the size-frequency 
distribution of dust, and the relative contributions of asteroidal and cometary material. The 
model can also be used to provide constraints on trace signals from other sources, such as 
dust associated with the ’’Plutinos” , objects captured in the 2:3 resonance with Neptune. 


Accomplishments 

• We now have a better appreciation for the relative importance of the sources that feed the 
cloud: as a result of our detailed models of the Solar System dust bands we believe that the 
cloud is predominantly asteroidal. Perhaps as much as 30% of zodiacal dust resides in the 
asteroidal dust bands alone (Figure 1), meaning that perhaps 90% of zodiacal dust originates 
in the asteroid belt and 10% comes from comets (Grogan et al. 2001, Dermott et al. 2003) 

• The dust band models also provide an estimate for the size-frequency distribution of dust 
particles in the main-belt, the result indicating that the dust band emission is dominated by 
large particles, greater than 100 microns in diameter (Figure 2, Grogan et al. 2001) 
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• We have shown how the orbital distribution of dust particles changes markedly over their 
dynamical lifetime, and that the action of secular resonances at the inner edge of the asteroid 
belt disperses low inclination, low eccentricity dust particles into the asteroidal background, 
making them appear more ’comet-like’ (Figure 3, Grogan et al. 2001; Kehoe et al. 2002a). 
This also has implications for the origin of extraterrestrial material accreted by Earth, and 
how it varies over time (Dermott et al. 2001; Kortenkamp et al. 2001) 

• We have shown how the stocastic nature of asteroidal collisions results in large-scale injections 
of collisional debris into the zodiacal cloud, with the result that the optical depth of the cloud 
can ’spike’ by orders of magnitude on numerous occasions over its lifetime (Figure 4, Grogan 
et al. 2001; Dermott et al. 2001) 

• We have provided constraints for the structure and detectabilty of dust particles trapped in 
a ’Plutino disk’, material captured in the 2:3 resonance with Neptune (Holmes et al. 2003) 

• We have further demonstrated how the planetary system imposes structure on the zodiacal 
cloud by warping its plane of symmetry, and by offsetting the cloud from the Sun (Dermott et 
al. 2001, Dermott et al. 2003). This information is directly applicable to the understanding 
of structure in Vega-type debris disks, and has led to the success of our subsequent NASA 
proposals (see Follow-on grants section). 


Follow-on grants 

Many of the results of this project, while directly applicable to the zodiacal cloud, are also relevant 
to the study of debris disks around other stars. A cornerstone of NASA’s Origins program is an 
understanding of the origin and evolution of planetary systems. Over the next decade observations 
from a variety of observing platforms, including SIRTF, Keck Interferometer, SIM, JWST, lead- 
ing to Terrestrial Planet Finder will allow us to detect extrasolar planets, and make meaningful 
estimates on how common these planetary systems might be. However, zodiacal emission in both 
our own system and in the target system may have the potential to seriously impact the chances 
for success of these missions. Addressing these concerns along with further studies of the zodiacal 
cloud are the subject of three follow-on grants, one already funded and two more which have been 
submitted this year. 

Funded 

• NASA Origins of Solar Systems, 2003-2005, ’’Dynamical Modeling and Detectability of Debris 
Disk Structure - Laying the Groundwork for Terrestrial Planet Finder” , PI Dr. Keith Grogan 

Submitted 

• NASA Terrestrial Planet Finder Foundation Science, 2004-2006, ’’The Nature of Asteroidal 
and Cometary Exo-zodiacal Clouds”, PI Dr. Keith Grogan 

• NASA Astrophysics Data Analysis Program, 2004-2006, ”MSX Constraints on the Origin and 
Structure of the Zodiacal Dust Bands”, PI Dr. Keith Grogan 
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Figure 1 : The 7 = 1.3 and 7 = 1.0 profiles in this figure were calculated assuming that the dust 
particles migrate from the asteroid belt to the inner solar system interior to 1 AU. The figure shows 
a comparison of the thermal emission obtained from two dust-band models to the corresponding 
unfiltered IRAS profile in the 25-/xm wave band and indicates that the dust bands contribute 
approximately 30% to the total near-ecliptic, thermal emission. Also shown (bottom profile) is 
the amplitude of the dust-band material calculated assuming that the dust-band material remains 
confined to the asteroid belt exterior to 2 AU. 
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Figure 2: Evidence for large particles dominating the dust band emission: dust band model residuals 
as a function of size-frequency index q. These residuals are calculated for a single line of sight 
denoted by the legend ’L, 304, 90’ meaning a longitude of Earth of 304° and a solar elongation of 90° 
in a direction leading the Earth in its orbit. They are found for the ten degree band by calculating 
the rms (observation-model) over two 5° wide latitude bins to cover the north and south bands for 
both the 12 and the 60 micron wavebands. The residuals are minimized around q = 1.4, although 
it becomes difficult to discriminate the exact value for low q since the particles in the model act 
like grey bodies (large diameter compared to the wavelength). The particle composition is taken 
to be astronomical silicate. 
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Figure 3: Vaxiation of the forced inclinations (representative of the gravitational perturbations of 
the planets) and the forced longitudes of ascending node with heliocentric distance for 10, 100 
and 200 micron diameter Eos family asteroidal dust particles at the present epoch. The dashed 
lines show the present osculating inclination and node for Jupiter, demonstrating that Jupiter 
dominates the orbital distribution of the particles beyond 2 AU. As the particles spiral towards 
the inner Solar System under the influence of Poynting- Robertson drag, a strong secular resonance 
near 2 AU strongly effects both the forced inclinations and nodes. 
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Figure 4: Variation with time of the total cross-sectional area of dust associated with the breakup 
of precursor asteroids that were big enough and numerous enough to supply all of the observed 
collision products of the main belt. This simulation modeled the stochastic breakup of asteroidal 
fragments, due to collisions, down to a diameter of 100 /im. A size-frequency distribution with a 
power-law index q — 1.9 was assumed for the initial breakup of each asteroid and this accounts for 
the height of the ’spikes’ in the plots 
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